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Abstract 

The evolution of a quantity, described by a function of space and time, relates 
the first derivative in time of this function to a spatial operator applied to the 
function. The initial value of the function at time t = is given. 

The fractional extension of this evolution consists of replacing the first deriva- 
tive in time by a fractional derivative of order a, < a < 1. 

We give a relationship between the solution of the equation of evolution and 
the solution of the equation belonging to its fractional extension. 
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1. Introduction 

Let u(x,t) be a function of space and time; x £ W 1 , t £ The equation, 
with initial condition, 

^l = [L(x)u](x,t) (1) 
u(x,0)=f(x) (2) 



is called an equation of evolution. The linear spatial operator L(x) depends on 
the problem at hand. We assume that the solution u(x,t) of this problem has 
been found. 
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The fractional extension of the above equation of evolution is given by the 
integral equation 



1 



u a (x, t) = u(x, 0) + — / dr(t- r)"" 1 [L(x)u a ] (x, r) (3) 



o 

< a < 1 . (4) 

For a = 1 we recover the solution of the equation of evolution. 

In this paper we find a relationship between the solutions u a (x, t) and u(x, t). 
A generalization of this relationship is also proposed. 

2. Evolution and its fractional extension I 

The equation of evolution Eqs. (H), (|2|) can be written as 

u(x,t) = f(x) + [ dr [L(x)u](x,t) (5) 
Jo 

u(x,0) = /(x). (6) 
The Laplace transform u(x,p) in the variable t is given by 

POO 

u(x,p) = / dte~ pt u(x,t) (7) 
Jo 

and leads to 

[L(x)u](x,p) =pu(x,p) - f(x) . (8) 
The fractional extension of the above equation of evolution reads 

1 f* 

u a (x, t) = f{x) + — - / dr(t- r) a - 1 [L(x)u a ] (x, t) . (9) 
r («) Jo 

The Laplace transform u a (x,p) in the variable t is given by 

/•oo 

u a (x,p) = / dte~ pt u a (x,t) (10) 
Jo 

and leads to 

[L(x)u a ] (x,p) = p a u Q (x,p) - f(x)p a - 1 . (11) 

Lemma 1. The solutions u(x,p) and u a (x,p) are related by 

u a (x,p) =p a - 1 u(x,p a ) . (12) 
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Proof. 



[L(x)u a ] (x,p) = p*- 1 [L(x)u] (x,p a ) (13) 
= p a - 1 {p a u{x,p a )-f{x)} (14) 
= p a u a (x,p)- f(x)p a - 1 (15) 



which is Eq. (11). 



We now look at the Mellin transforms of u(x, t) and u a (x, t) in the variable t. 

poo 

u(x,s) = / dtt s ~ l u{x,t) (16) 

J 

u a (x,s)= / dtt s ~ 1 u a (x,t) . (17) 
Jo 

The relationship between the Mellin and Laplace transform of a function cp(t) 

is 

i r°° 

' dpp- s 4>(p) . (18) 



o 



" v ' r(i-s 

Lemma 2. u(x,s) and u a (x, s) are related by 



&*{x, s) = - r p ^ a } u (x, -) . (19) 
a 1 (1 — s) \ a) 



Proof. From Eqs. (|TsD , ([[D we get 

U a {x, S ) = — r / dpp- S p a - l u{ X ,p a ) . (20) 

r (! - s ) Jo 

The change of variable 

P a = q (21) 

leads to 

11 f°° _ ± 

u a (x,s) = -— -/ dqq «u{x,q) (22) 

aT(l - s) J 

or 

U a (x, s) = — ~^7T~ T r f 1 - -) « ( x i -) • ( 23 ) 

a 1 (1 — s) \ a) \ a) 

This is Eq. (U). 
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3. Generalized Mittag-Lefner Functions 



The generalized Mittag-Lefner functions pi] are given by 



oo 

F a ,(,)=r(/?)p-i) fc — 



(/? + ak) 



with 



z>0,0<a<l,/3>a. 

= 1 give the usual Mittag-Leffler functions which we denote by F a (z). 
The H- function representation of F a p{z) j|, ||, Q is given by 



(0,1) 
(0,l)(l-/3,a) 



(24) 



(25) 



(26) 



The generalized Mittag-Leffler functions are the Laplace transform of probability 
measures on The corresponding probability densisties are given by 



f aP { Z ) = nmii z 



((3 — a, a) 
(0,1) 



(27) 



(3 = 1 gives the probability density of the usual Mittag-Leffler functions which is 
denoted by f a (z). It has the H- function representation 



f a {z) =H\\[z 



1 — a, a 
^ (0,1) ' 



< a < 1 



(28) 



and 



h(z)=5(z-l) 



(29) 



For < a < 1, f a (z) is an entire function and vanishes exponentially for large 
positive It has the power series representation 



oo _. 



fc=0 

< a < 1 , zeS + . 



(30) 
(31) 



4. Evolution and its Fractional Extension II 
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The Mellin transform of the probability density f a (z) is given through Eq. (28) 



as 



fa{s) 



T(s) 



(32) 



by 



r(l — a + as) 

Lemma 3. The inverse Mellin transforms of u(x,s) and u a (x,s) are related 

(33) 



u a (x,t) = / dz f a (z)u (x , t a z) 
Jo 



Proof. Let M denote the performance of the Mellin transform. From 
we have 

r»0O 



M \ dzf(tz)g(z) 



o 



(s) = f(s)g(l-s) 



[Mf(at a )](s) = -a-if(-) 



M I dzf(z)g(t a z) 
o 



(s) = - 



M / dzf{z)g{tz) 



From Eq. (|32| ) we find 



and thus from Eq. ( |19| 



This leads to 



M 1 a) ~ r(i-T) 



u a (x, s) = -f a ( 1 - -) u (x, -) . 

a V aJ \ aJ 



1 



u n (x, s 



a 



M / dz f a (z)u(x,tz) 
Jo 

/•oo 

M / dz f a (z)u(x,t a z) 
Jo 



to 



Thus 



U a (x,t) = / dz f a (z)u(x,t a z) 



(34) 
(35) 
(36) 

(37) 

(38) 

(39) 
(40) 

(41) 



Theorem 1. The solution u(x, t) of the equation of evolution Eq. ([|) and the 
solution u a (x,t) of its fractional extension Eq. (|J) are related by 

u a (x,t) = t~ a dzf a (t- a z)u(x,z). (42) 
Jo 



Proof. Make the variable transform t a z = w in Eq. (|33|). 
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5. The Green's functions 

The solution of the equation of evolution Eq. (||) can be written as 

u(x,t) = J dyG(x,y;t)f(y) (43) 

where we assume that the Green's function G(x, y; t) is known. 

The solution of the equation of its fractional extension Eq. (|9|) can be written 

as 



u a (x,t) = J dyG a (x,y;t)f(y) . (44) 
According to Theorem 1, Eq. (|42|) the two Green's functions are related by 

POO 

G a (x,y;t)= t~ a / dz f a (t~ a z) G(x, y; t) . (45) 
Jo 

6. Applications 

1. The fractional diffusion equation ||. The equation of evolution is given by 

du(x, t) 



dt 

where 



Au(x,t) (46) 



n ^2 

A = ^I? (47) 

k=l ° X k 

and the initial condition 

u(x,0) = f(x). (48) 

The Green's function for this problem is 

G(x,y;t) = G(\x-y\,t) (49) 

where with 

r = \x — y\ (50) 
G(r,t) = (4vrt)~f e -^ . (51) 
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The Green's function for the fractional extension of Eq. ( f46| ) is then 

/>oo 

G a (r, t) = t~ a / dz f a {t~ a z) G(r, z) . (52) 
J o 

We compute this integral by looking at its Mellin transform in the variable 
t. 

<5.(r,.)-i/.(l-i)«S(r.i) (53) 
- ir ^-i)6(r,l). (54) 



a r(l — s ) \ a 

From 1| we get 



and thus 



G(r,s) =7r-2 2-^r is - n r(--s) (55) 



a V 2 aJ 1 (1 — s 



- 2 ' iX^' 2^ 



(57) 



Inverting the Mellin transform leads to 

G a M)=7r-f2-V-^?° Qrt-f 

This is indeed the Green's function as given in 

2. Fractional Black-Scholes equation Q 

The equation of evolution (Black-Scholes equation) is given by 

S > , t < T (59) 

and boundary condition 

C(S,T) =max{S - E,0} . (60) 
With the transformation 

t = T — %t , C(S, t) = A(S, r) (61) 

(J* 

2r 

Ao = — (62) 
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we have the equation of evolution 



with the initial value 



^_^A BA XA 



A(S,0) =max{S - E,0} 



The solution is given by 

A(S,t) = SN(dx) - Ee- XoT N(d 2 ) 

where 



d x = (2r)-3 
d 2 = (2r)-5 



and 



N(d) 



In ( | ) + (A + l)r 
ln(|) +(A -l)r 



dx e 2 



(63) 

(64) 

(65) 

(66) 
(67) 

(68) 



According to Theorem 1 Eq. ( |4^ ) the solution of the fractional extension of 
the Black-Scholes equation 



1 



A a (S, r) = A(S, 0) + — - / dz(r-z) 



\a-l 



r(a) 

2 d 2 A,(S,z) , x C £U Q (S» 



s 



dS 2 



+ \ S 



dS 



X A a (S, z) 



is given by 



roo 

MS, r) = r- a / dz f a {r' a z) A(S, z) 
Jo 



(69) 



(70) 



7. Generalization 

According to Eq. (^) the solution of the equation for the fractional extension 
an equation of evolution is given by Eq. (|33|) 

POO 

u a (x,t) = / dz f a (z)u(x,t a z) (71) 
Jo 
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f a (z) is a measure on 5? + and fi(z) = 5(z — 1). 

The Laplace transform of f a is the usual Mittag-Lemer function. 

The generalized Mittag-Lemer function Eq. (24) is the Laplace transform of 
the measure f a /3(z) on Eq. (27). We have 



f a/3 (x) = T((3)Y / (-l) 



k=0 



T{(3 — a — ak) k\ 



a < 1 , (3 > a 



(72) 



{(3 - 1)(1 - xf~ 2 , 0<x<l 







1 < x 



, p > 1 (73) 
fn(x) = S(x - 1) . (74) 
It is thus tempting to generalize from the solution u a (x,t) in Eq. (fnl) to 



u a/3 (x,t) 



dz f af i{z)u a p(x,t a z) . 



(75) 



This expression uses the form factor f a p and will be studied elsewhere. Other 
form factors might also be considered. 



8. Summary 

We assume that we know the solution of a general equation of evolution. The 
solution of the equation of its fractional extension is then related to the above 
solution by what we call a form factor. This form factor belongs to the Mittag- 
Lemer function. A fractional generalization is proposed by using the form factor 
belonging to the generalized Mittag-Lemer function. 
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